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We report on the existence of nonlinear surface waves which, on the one hand, do not require 
the threshold energy flow for their excitation, and, on the other hand, extend into media at both 
sides of the interface at low powers, i.e. can not be reduced to the conventional Tamm states. Such 
waves can be excited if the refractive index in at least one of the materials forming the interface 
is periodically modulated, with properly selected modulation depth and frequency. Thresholdless 
surface solitons can be stable in the entire existence domain. 

PACS numbers: 



The history of the surface waves goes back to the sem- 
inal paper of Tamm who upon consideration of the 
Kronig-Penney model for an electron in a crystal, has 
shown that at the interface between a periodic and homo- 
geneous structures there can exist localized states, which 
cannot be supported by a surface between two homo- 
geneous media. Four decades later, in Ref. Q Tamm's 
ideas were introduced in the optics of periodic structures, 
where localized states of the electric field at the bound- 
ary between a homogeneous and stratified media were 
predicted. At the same time it was understood Q that 
the nonlinearity can dramatically change the situation 
leading to new surface states, whose analogs do not ex- 
ist in the linear structures. Nowadays the issue of the 
nonlinear surface modes attracts growing attention (see 
e.g. the recent review [3| and the references there in). In 
particular, it was suggested that nonlinear surface modes 
may form at the interface of periodic and uniform non- 
linear media Such states were thoroughly studied 
and observed in both focusing and defocusing 0-Q 
media. 

One of the main features of the nonlinear surface states 
at the boundary between two semi-infinite media, stud- 
ied so far, stems form the fact that they have no linear 
limit: such modes possess the energy flow threshold be- 
low which they do not exist. Such thresholds have been 
widely discussed in the literature, see e.g. So far, 

thresholdless surface states were observed only at the in- 
terface of two periodic structures separated by a suitable 
interval El , or at the interfaces of lattices with de- 
fect surface channel [IH, [TJ] , or at the interfaces of non- 
periodic lattices [l4j. Notice however, that in all these 
settings surface modes bifurcate from well-localized lin- 
ear defect states (i.e. from the Tamm states). 

In the present Letter we show that there exists one 
more type of the nonlinear surface modes, which on the 
one hand do not require threshold energy flow for their 
excitation and on the other hand do not have linear lo- 
calized counterparts, unlike all surface modes reported so 
far. 

We consider beam propagation which can be described 



by the equation for the dimensionless electric field q: 

iqe, = -%n + R{v)Q + (J \ ( l\ 2( l- (!) 

Here £ and r] are the longitudinal and transverse coordi- 
nates normalized to the diffraction length and the char- 
acteristic beam width. It will be assumed that the plane 
i] = separates two media. The right hand side medium 
(rj > 0) is considered periodic with the modulation of 
the refractive index described by R r i g ht(v) = Rr cos(2?y) 
(we emphasize that the particular choice of the cos-like 
refractive index is for the illustration purposes only, and 
the results reported below are valid for any other type 
of refractive index modulation). For the refractive in- 
dex of media on the left hand side, 77 < 0, we sup- 
pose Ri e ft(rj) = Rifi + cos(2?y) and consider the two 
different examples: of the homogeneous Kerr medium 
(Rl,i = 0) and periodically modulated one, Ri t i ^ 0. 
We concentrate on the cases of two focusing, a = —1, or 
defocusing, a = 1 media (the generalization of our argu- 
ments to the case of the interface between focusing and 
defocusing media is straightforward). 

To explain the main idea we start with the case, where 
Rieftiv) = Ri,o- We are interested in stationary solutions 
q = e lb ^w(ri) (here b is the real propagation constant and 
w(r]) describes the soliton profile) bifurcating from the 
linear spectrum. In order to get insight into properties 
of such states it is instructive to start from the linear limit 
q — » (or formally a = 0). To this end we concentrate 
on the linear equation ~bw(rj) — —w m + R(rj)w(r)), ob- 
tained from JT]) by neglecting the nonlinear term, which 
at 77 — > 00 decays more rapidly than the linear ones 
(notice that 10(77) — > 10(77) at 77 — > 00). Now we re- 
call the standard arguments giving the Tamm mode (see 
e.g. [Hj]). The linear solution for 77 < is given by 
wi(rf) = C/exp(/^/77) where /// = ^/i?;,o + b, while the 
solution at rj > has the form of the Bloch function 
w r = P(rj) cxp (—fj, r rj) where P(j]) is a it- or 27r-periodic 
function and fi r is the Floquet exponent (we have taken 
into account that the propagation constant b belongs to 
one of the gaps of the spectrum of the lattice Rrightiv))- 
Then the condition for existence of the Tamm mode in 
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FIG. 1: (Color online) Schematic representation of interfaces 
between homogeneous and periodic media (left panel) and 
between two different periodic media (right panel). Shad- 
owed areas show the linear spectrum (allowed bands). The 
bold parts of the axis b indicate the intervals to which the 
propagation constant must belong for existence of any surface 
soliton. The bold horizontal line shows the edge of the lin- 
ear spectrum, from which the thresholdless solitons bifurcate. 
The dashed lines show examples of location of the propagation 
constants of the surface solitons. The values of b indicated in 
the left and right panels correspond to the specific examples 
shown below in Figs. [5] and 2] 



the general case reads /i r = P v (0)/P(0) — m. In our par- 
ticular case, where max, R r ight(jl) = Rright{0), we have 
P v (0) = and thus /Lt r — > 0, at b — > —Ri.o- In this limit 
we also have that w r , v (0) — > (hereafter the subscript T) 
stands for the derivative, i.g. wi tV = dwi/drf). Since, in 
a periodic medium, the Floquet exponent tends to zero 
only if the propagation constant approaches one of the 
gap edges, we conclude that the limit fi r — > is possible 
only if Rrightiv) is such that —Rifi coincides with one 
of the boundaries of the gaps (or allowed zones). This 
leads us to the first condition, necessary for existence of 
the thresholdless surface solitons: the boundary of the 
6-domain where the localized states exist in the homoge- 
neous medium must coincide with one of the boundaries 
of the stop gaps of the periodic medium. This situation 
is illustrated in Fig.QJa). (If P^O) ^ 0, i.e. the interface 
between the two structures does not coincide with one 
of the extrema of the periodically modulated refractive 
index, the respective corrections must be introduced in 
the consideration.) 

Turning now to the nonlinear case, we consider the 
limit b — > —Rifi (illustrated by the dashed lines in Fig.[T]) 
when surface solitons are very wide and their envelopes 
inside right (periodic) medium can be rather accurately 
described by the homogeneous nonlinear Schrodinger 
(NLS) equation for the envelope of the Bloch wave, with 
the effective diffraction D = — d 2 b n (k) /dk 2 \k=o.i where 
b n (k) is the n-th band of dispersion relation of Rrightiv) 
(k is the Bloch wavevector; k = and k = 1 correspond 
to the center and the boundary of the Brillouin zone). 

First we consider the case of the defocusing media 
(cj = 1). Then soliton can exist only if D < and 
its maximum can be located only in the right (peri- 



odic) medium. For 77 < the solution is given by 
ui[ = v2"mz/ sinh (/x/ (77; — 7/)), where rji > is a constant, 
which must be found from the continuity of w and w n at 
the surface. When 77 = 0, the function wi(rj) (defined at 
77 < 0) achieves its maximum w^ max = V^fii/ sinh (fJ-iTii). 
In this point one verifies wi jV <~ fj, 2 / sinh 2 (/z; rji) ~ wf 
for fi r — > 0. On the other hand, for 77 > 0, where 
D > 0, the envelope of the Bloch state is described by 
the cosh _1 -like function, and at 77 = one has the re- 
lation w r< r,(0) ~ w r (0) for /i/ ~ — > 0. Since, there 
must also verify u>;(0) = w r (0) and wi tV (0) = w r:V (0), 
the above asymptotics for wi jr are not compatible. This 
means that the amplitude of the mode does not go zero, 
i.e. no thresholdless surface solitons can exist. 

Consider now the case of the focusing media (a = 
— 1). Now for 77 < the solution is given by wi = 
\Pl\xiJ cosh (/i/(ry — 77/)). For 77 > one must require 
now D < (notice that, if D > 0, the solution in the 
periodic media is of sinh^Mike type, and the above ar- 
guments about existence of thresholdless solitons at the 
interface between focusing and effectively defocusing me- 
dia, inhibit the existence of thresholdless modes). Now 
the asymptotics of the left-medium and right-medium so- 
lutions are compatible, i.e. the limit b — > —Rifi can al- 
low for existence of the solution localized about the sur- 
face. Thus we conclude, that the solution, localized in the 
vicinity of surface, can exist when D > 0, what is satisfied 
for all upper (lower) boundaries of the allowed (forbid- 
den) zones of the lattice. With this assertion we now 
can formulate the requirement for the existence of the 
thresholdless surface solitons in a focusing medium: the 
boundary of the existence domain of the localized states 
in the homogeneous medium must coincide with one of 
the upper boundaries of the allowed bands. 

The numerical analysis of this case is illustrated in 
Fig.[5Ja), where the branch of the solution is depicted in 
the plane (b, U) where U — J_ \w\ 2 drj is the total energy 
flow along the surface. We consider the case correspond- 
ing to Fig. 1(a) when the lower boundary of existence 
domain inside the uniform medium b = —Ri t o coincides 
with lower boundary of semi- infinite gap in right medium. 
One observes that the soliton branch indeed bifurcates 
from the continuous spectrum, i.e. surface solitons can 
be excited at any infinitely small power. The surface 
soliton width decreases with the increase of soliton en- 
ergy flow, as it follows from the comparison of Figs. &h) 
and[2Jc). Such solitons are completely stable (we observe 
that for them dll/db > 0, similar to the stability of soli- 
tons in a homogeneous medium ensured by the Vakhitov- 
Kolokolov criterion). The same situation occurs when 
propagation constant falls into first gap of periodic right 
medium, when the — i?/,o coincides with the upper edge 
of the second allowed band (see Fig(3J) . These solutions 
however are unstable. 

Let us now turn to a boundary between two dissimilar 
lattices, which we studied for the focusing case a = — 1. 
Employing the above arguments one can conclude that 
for existence of thresholdless surface solitons the coinci- 
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FIG. 2: (Color online) (a) Energy flow vs propagation con- 
stant for the interface between periodic R r = 3 and homoge- 
neous -R; : o = —0.9368 media for the focusing (a — —1) non- 
linearity. Shaded area corresponds to the first allowed band 
of the lattice spectrum; (b,c) shape of surface soliton at cor- 
responding points of panel (a) for b = 0.94 (b) and b — 1.7132 
(c). The zone structure for situation is schematically depicted 
in the left panel of Fig. [T] 



pie of soliton family whose propagation constant emerges 
from semi-infinite gaps (with coinciding lower edges) of 
both lattices is shown in Fig. |4] Such solitons are stable 
too in the entire domain of their existence. 

To conclude, we have reported a possibility of existence 
of the thresholdless surface solitons, which in the linear 
(i.e. small energy flow) limit are transformed into small- 
amplitude extended states approaching Bloch waves in- 
side the periodically modulated medium and having long 
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FIG. 3: (Color online) (a) Energy flow vs propagation con- 
stant for the interface between periodic R r — 3 and homoge- 
neous Rifi = 2.166 media for the focusing nonlinearity. Left 
and right shaded areas correspond to the first and second 
allowed bands; (b) and (c) shapes of surface solitons at the 
corresponding points of panel (a) for b = —2.15 and 6 = 0. 

dence of two upper boundaries of the allowed bands of the 
left and right lattices is required. A representative exam- 



FIG. 4: (Color online) (a) Energy flow vs propagation con- 
stant for the interface between periodic media with parame- 
ters R r — 3, Rlo = —0.815, Ria = 1 media for the focusing 
media. Shaded area corresponds to the first band of the spec- 
trum; (b,c) shapes of surface solitons at corresponding points 
of panel (a) for b = 0.9456 (b) and b = 1.7183 (c).The zone 
structure for this situation is schematically depicted in the 
right panel of Fig. [T] 



slowly decaying tails inside the uniform medium. In this 
sense the reported surface modes differ drastically from 
conventional surface waves possessing excitation thresh- 
olds and from modes reducing in the linear limit to lo- 
calized Tamm states. 

The main condition for existence of the thresholdless 
solitons is matching of the boundaries of the linear spec- 
tra at the both sides of the interface. In practical terms 
the required conditions can be implemented in various 
ways, e.g., by adjusting of geometrical properties of the 
lattices (periods, depths of modulations, etc.), by using 
temperature gradients (i.e. change of the refractive index 
by heating one of the two media), by doping of the ma- 
terials that results in modifications in refractive index, 
etc. 
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